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We provide an effective description of fractional topological insulators that include the fractional 
quantum spin Hall effect by considering the time-reversal symmetric pendant to the topological 
quantum field theories that encode the Abelian fractional quantum Hall liquids. We explain the 
hierarchical construction of such a theory and establish for it a bulk-edge correspondence by deriving 
the equivalent edge theory for chiral bosonic fields. Further, we compute the Fermi-Bose correlation 
functions of the edge theory and provide representative ground state wave functions for systems 
described by the bulk theory. 



I. INTRODUCTION 

Laughlin initiated the theoretical exploration of the 
fractional quantum Hall effect (FQHE) by proposing 
wave functions for the ground states of interacting elec- 
trons in the lowest Landau level at filling fractions 
v = l/(2m -1-1), m g Z.El The experimental observa- 
tion of a plethora of fractional Hall plateaus at other 
filling fractions lead to the construction of a hierarchy 
of wave functions out of Laughlin's wave functioUjI^HIl 
and the development of the composite fermion picture.'^ 
These approaches were later reconciled, and unified by 
the effective description of the FQHE in terms of multi- 
component Chern-Simons theories in (2 -I- l)-dimensional 
space and time.'SESl These topological effective theories 
for the hierarchy of the FQHE deliver a correspondence 
between the physics in the two-dimensional bulk and the 
physics along one-dimensional boundaries at which the 
two-dimensional sample terminates.'^SHli] 

It is possible to double the Chern-Simons effective the- 
ory representing the universal properties of the FQHE at 
some filling fraction v — l/(2m -1-1), m S Z so as to ob- 
tain a time-reversal symmetric theory. This approach has 
been used to interpret a fully gaped superconductor as 
an example of a topological phase,'^^!!!!] g^j^jj „ more gen- 
erally - to explore the universal properties of interacting 
theories with an emergent local gauge symmetry (see 
Refs. that si gna ls the phenomenon of spin and 

charge separation.'SHEII 

A more urgent impetus for the construction of effec- 
tive time-reversal symmetric topological field theories in 
(2 -I- l)-dimensional space and time arose with the theo- 
retical prediction of time-reversal symmetric topological 
band insulators, shortly followed by their experimental 
discovery.'S^'l^ll These band insulators realize the counter- 
parts to the integer quantum Hall effect and their discov- 
ery suggests that a time-reversal symmetric counterpart 
to the FQHE might emerge from interacting itinerant 
electrons in a crystalline environment. 

From the outset, this endeavor follows a different line 
of logic than the FQHE, as it is not based on pre-existing 



experimental evidence. Past experience with the FQHE 
has thus guided recent attempts to either construct time- 
reversal symmetric edge theories or to construct time- 
reversal symmetric bulk wave functions supporti ng loc al 
excitations carrying fractional quantum numbers.'^'^^Hill 
While numerical support for a time-reversal symmet- 
ric topological phase of matter was given by Neupert 
et al. in their study of a lattice model for interacting 
itinerant electrons,H2 a description in terms of an effec- 
tive theory is desirable to reveal the universal proper- 
ties of such a phase. In Ref. SH the universal proper- 
ties such as the topological degeneracies of the ground 
state manifold were explored with the help of a family of 
edge theories. In this paper, we are going to construct 
the corresponding bulk topological theory by generaliz- 
ing the hierarchy of Abelian FQHEs to the hierarchy of 
Abelian fractional quantum spin Hall effects (FQSHEs) 



in Sec. Hll We will show in Sec. HI the correspondence 



between the bulk theory and the edge theory whose sta- 
bility to the breaking of translation invariance and resid- 
ual spin-1/2 U(l) symmetry was studied in Ref. Fi- 
nally, we shall generalize in Sec. IV the wave functions 
supporting the Abelian FQHE for a fractional filling of 
the lowest Landau level to wave functions supporting 
an Abelian FQSHE. These time-reversal symmetric wave 
functions are built from the holomorphic and antiholo- 
morphic single-particle wave functions belonging to the 
lowest Landau level when the applied uniform magnetic 
field is pointing down or up, respectively. For the reader 
who wants to skip the derivations, we provide a detailed 
summary of our results in Sec. |V] 



II. TIME-REVERSAL SYMMETRIC ABELIAN 
CHERN-SIMONS QUANTUM FIELD THEORY 

Let us start by summarizing some of the results that 
we will derive in this section. We shall construct a class of 
incompressible liquids, each of which is the ground state 
of a time-reversal symmetric (2+l)-dimensional Chern- 
Simons quantum field theory that depends on 2A^ flavors 
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of gauge fields ^{t, x), where i = 1, . . . , 2N labels the that leaves the action (2.2b) invariant. In this represen 



flavors and /i = 0, 1, 2 labels the space-time coordinates 
= {t,x), with the action 



S := 



c s 

Qt \ 9^ 



27r 



2^ ^t^^u^Pj 

(2.1a) 



Here, K^j are elements of the symmetric and invertible 
2N X 2N integer matrix K. The integer-valued com- 
ponent Qj of the 2iV-dimensional vector Q represents 
the z-th electric charge in units of the electronic charge 
e, which couples to the electromagnetic gauge poten- 
tial a;). Similarly, 5*^ is an integer- valued com- 
ponent of the 2iV-dimensional vector S that represents 
the i-th spin charge in units of s associated to the up 
or down spin projection along a spin-1/2 quantization 
axis, which couples to the Abelian (spin) gauge potential 
B^{t, x). The operation of time reversal maps a- ^{t, x) 
into —g'^'' a-^_^^pf ^{—t,x) for i = 1, • • • , iV and vice versa. 
Here, g^^, := diag(+, — , — ) = g'^'^ is the Lorentz met- 
ric. In Eq. (2.1a) x e fi, where C is a region of 



two-dimensional Euclidean space, which for the discus- 
sion of the bulk theory we consider to have no boundary, 
dfl — 0. The domain of integration M is unbounded in 
time t. We will show that time-reversal symmetry im- 
poses that the matrix K and the vectors Q and S are of 
the block form 



K = 



K A 

AT -K 



Q = 



S = 



Q 



(2.1b) 



with Q an integer iV-vector, while n — kJ^ and A ~ — A^ 
are symmetric and antisymmetric integer-valued N x N 
matrices, respectively. 

The doubled structure of the theory is even more evi- 
dent if we express it as a BE theoryjE^f^i.e., by defining 



1 



i = l,...,iV, (2.2a) 
for /i = 0, 1,2. This basis allows to re-express the effec- 



tive action (2.1a 



S := / dtd^x 



1 



(-) 



+ A^d,a\+> + - p; B^d.al; 

(2.2b) 

In this representation, the indices in sans serif fonts i,j 
run from 1 to N. The coupling between the pair of gauge 
fields a'^'^^ and a^~^ is off-diagonal in the BE labels ±. 
This is a consequence of time-reversal symmetry, which 
is implemented by 



it,x)^Tg''''al^>i-t,x) 



(2.2c) 



tation, the electromagnetic gauge potential A couples to 
the +-species only, while the spin gauge potential B cou- 
ples to the — species only. The N x N integer-valued 
matrix }i in the BE representation is related to the block 
matrices k and A contained in K from Eq. (2.1b) through 



= K - A. 



(2. 2d) 



The degeneracy of the ground state is obtained for ei- 
ther description, i.e., the one in terms of the flavors 

with i — 1, - ■ ■ , 2N or the one in terms of the flavors ap-* 
with i = 1, • • • , TV, from 



GS 



det 



X 
>/ 



(det >f)^ 



(2.3) 



If the underlying microscopic theory describes fermions 
with a residual spin-1/2 U(l) (easy plane XY) symmetry, 
it is then meaningful to define the quantized spin Hall 
resistance 



2tt 



X 



(2.4a) 



The filling fraction is here defined so that it is unity for 
the integer quantum spin Hall effect and therefore given 
by 



= Q^ X^^ Q. 



S 



(2.4b) 



We now turn to the hierarchical construction of the 
states described by this quantum field theory. As a warm- 
up, we begin by reviewing how a one-component Chern- 
Simons quantum field theory in (2-|-l)-dimensional space 
and time is related to the quantum Hall effect. We then 
construct recursively the multi-component Chern-Simons 
quantum field theory in such a way that it respects time- 
reversal symmetry. 



A. Brief review of the one-component 
Chern-Simons theory 



We start from the Lagrangian density 
in (2 + l)-dimensional space and time with the action 



5cs / / d^xCcs (2.5b) 



and partition function 

Zcs[A] j PHe^-^cs. (2.5c) 
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The dimensionless integer p is positive. The electro- 
magnetic coupUng (electric charge) e is dimensionfuU. 
It measures the strength of the interaction between an 
external electromagnetic gauge field A with the compo- 
nents = (yl", A) and a dynamical gauge field a with 
the components = {a'^,a). The symbol I? [a] repre- 
sents the measure of all gauge orbits stemming from the 
Abelian group U(l). 

The operation T for reversal of time is defined by 



A(t,x)^ +g^^'^A^{-t,x), 



(2.6a) 



(2.6b) 



for /X = 0, 1,2. We also posit that T is an anti-unitary 
linear transformation. If so, one verifies that C^g is odd 
under reversal of time. 

Define the electromagnetic current to be the 3-vector 



J; 



^ _ 1 SS^g 



cs 



h SA„ 



(2.7a) 



2nh 



for ^ = 0, 1,2. Because the Levi-Civita tensor with the 



component e 
is conserved, 



012 



1 is fully antisymmetric, this current 



Now, the equations of motions 



(2.7b) 



= '^'^cs 
5a,, 



i^e^''^d„a, + ^e^^''^d,,A, 



2n^ -■^^'"2^^- "''^^ (2.8) 



can be used in conjunction with Eq. (2.71 to yield the 
conserved electromagnetic current 



cs 



if! 

p h 



ef""^ d„ A 



(2.9) 



and partition function 



Zl^\A,B\ := 1 V\a 



(+) „(-)! 



(2.10c) 



Equation ( 2.10 1 is a BE theory made of two copies of the 



Chern-Simons theory (2.5 1 with the specificity that the 
integer p enters with opposite signs in the two copies. We 
have also introduced two external gauge fields A and B 
with the couplings e and s, respectively. Eor the gauge 
field A, e will be interpreted as a total U(l) charge. For 
the gauge field B, s will be interpreted as a relative U(l) 
charge. If the underlying microscopic model is built from 
itinerant electrons, the gauge field A is the U(l) electro- 
magnetic gauge field that couples to the conserved elec- 
tric charge whereas the gauge field B is the U(l) gauge 
field that couples to the conserved projection along some 
quantization axis of the electronic spin, i.e., s = 1/2. 

This theory is invariant under the operation of time 
reversal defined by the anti-linear extension of 

af\t,x)^ T.g'^''a(±)(-t,a;)^Ta^±)^(f,i), (2.11a) 



Ait,x) ^ +g^'•' A^{-t,x) = +A^(i,i), (2.11b) 



B^{t,x) ^ -g^^'' B^{-t,x) = -i?''(i,A), (2.11c) 

for /X — 0,1,2. The component A^ of the external elec- 
tromagnetic gauge field A is unchanged whereas its vec- 
tor component A is reversed under reversal of time, just 
as the vector components of a^~\ This behavior is re- 
versed for the components of the external gauge field B 
that couples to the conserved U(l) spin current and the 
gauge field a'^+^. 

Since this theory is equivalent to two independent 
copies of the Chern-Simons theory (2.51, there are two 



independent conserved currents of the form (^2.7L 



which allows us to identify the filling fraction v ^ p ^ 
in this simple example, so that the quantum Hall con- 



ductance is given by cth 
units in which h= 1. 



From now on, we adopt 



(2.12) 



for /i = 0, 1, 2. Their transformation laws under reversal 
of time are 



B. One-component BF theory 

We start from the Lagrangian density in (2 + 1)- 
dimensional space and time 



TT 



with the action 



-BF 



(2.10a) 



(2.10b) 



(2.13) 



for p = 0,1,2. If the microscopic model is made of itin- 
erant electrons, we can thus interpret J!^ as the charge 
current and, if the model has a residual U(l) rotation 
symmetry of the electronic spin, represents the con- 
served spin current. The equations of motions 



= ^'^BF^' 



5a 



(±) 



(2.14a) 



for the dynamical compact gauge fields ' and a^^\ 
respectively, deliver the relations 



(2.14b) 
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and 



P 



(2.14c) 



for /i = 0,1,2, respectively. We conclude that, on the 
one hand, the charge current obeys the Hall response 

= 2s X —e^"'^ B., (2.15a) 

with fj, — 0,1,2 while, on the other hand, the spin current 
obeys the Hall response 



2e X 



27rp 



(2.15b) 



with = 0,1,2. 



for i = 1, • • • , n is also assumed. 

The operation of time reversal is the rule 

x^^i'^:=-g^,x-' (2.18a) 
together with the anti-linear extension of the rules 

al^^''{x)^Tg,^al^^-{£), (2.18b) 

for fi = 0,1,2 and i = 1, • • • , n that leaves the Lagrangian 
density (2.16c) invariant. 



The level n + 1 of the hierarchical construction posits 
the existence of the pair of quasiparticle 3-currents j^. 
that are conserved, i.e., 



(2.19) 



C. Time-reversal symmetric hierarchy 

The generic structure of the hierarchical construction 
is the following. Let n > be any positive integer. Define 
at the level n of the hierarchy the quantum field theory 
with the partition function 



Zj^^[A,B] := jv 



where the action is 

R n 

and the Lagrangian density is 

n 



i=i 



,a. 



(-) 

n 

(2.16a) 
(2.16b) 



+ ^^^H,M.A^^5^„(+) (2.16c) 



It also posits the existence of some even integer p„^i 

and 2n integers ll with i = 1, • • • ,n such that the 
constraints 



■?±,ri+l 



np, 



n+l 



(±) 

A 



(2.20) 



for /i — 0,1,2 hold. The constraint (2.201 means that 



any pair of flux quanta, arising when af^'* and ^ each 
support a vortex, creates a quasi-particle with charge 



2 ;[^VPn+i and spin 2 ll V-P«+i for i = 1, • • • ,n. 

This construction can be achieved from the partition 
function 



(-), 



,(+) 



'"n+l'"l ! 



X e ""+1 , 



with the action 



'■n+l 



(2.21a) 



5™f := jdtj d^xC^+, (2.21b) 



(-) 

A ■ 



Here, the dynamical gauge fields a^^^ are the n-tuplet 
with the components 



(2.17a) 



Moreover, the nx n matrix >f(") is invertible and has, by 
assumption, integer-valued matrix elements. The charge 
vector £1*^"-' has the integer-valued components 



= (1,0,--- ,0^ eZ" 

Finally, the compatibility condition 



(2.17b) 
(2.17c) 



and Lagrangian density 

/-TRS /.TRS 



_ Pn+1 ui^A (+) a (-) 



+ 



l^p.Ay,(-)^(-)g^^^(+) 

TT ^ ^ 



^,11, '^i^"n+l,A- 



(2.21c) 



Indeed, we can then define the conserved quasiparticle 
currents of type n to be 



,M 1 MfA Q (±) 



(2.22) 
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for /i = 0, 1, 2 and use the equations of motion 







da, 



(T) 



Pn+l uuXc, (±) _ 
TT ^ ' TT 



(2.23) 



9 a<^^ 



obeyed by the dynamical gauge fields a-^+i^i estab- 
lish that they indeed obey the constraints imposed in 
Eq. ((2:201 . 



Observe that if we introduce the two (n + l)-tuplets 



7(±) 



given by 



(2.24a) 



of dynamical gauge fields, then the Lagrangian £™i de- 
fined in Eq. (|2.21c ) takes the same form as >C™^ defined 
in Eq. ( 2.16c[ ) after the substitution n — > n -|- 1. The 
(n -|- 1) X (n -|- 1) matrix is then given by 



An) 



(2.24b) 



The {n + l)-component charge vector is given by 



(n+l) 



(f,0,-- - ,0)"^ e Z"+\ 



(2.24c) 



thus imposing a vanishing coupling of the external gauge 
fields A and B to a^j^\. The compatibility condition 



i-y 



(2.24d) 



for i = 1. • • • , n + 1 holds if and only if the integer p^-^^i 



The representation (2.24) is called the hierarchical rep- 
resentation. 



The operation of time reversal obtained from Eq. (2.18 ) 



by allowing i to run from 1 up to n -I- 1 leaves the La- 
grangian of level n -|- 1 invariant. Therefore, we have 
constructed a hierarchical time-reversal symmetric BE 
theory. 



D. Equivalent representations 



We define an equivalence class on all the actions of the 



form (2.2b I when there exists a linear transformation W 



with integer valued coefficients and unit determinant 
such that 



and 



x=W^ x'W 



g = W^ g', 



(2.25a) 



(2.25b) 



between any two given pairs (>f, g) and (>r', g') within an 
equivalence class. 



Example 1: The lower-triangular transformation 



(1 ... 0\ 

1 -1 ... 



V 



(2.26a) 



-V 



relates the hierarchical basis characterized by the charge 
vectors 



e = (1,0,... ,0)T 



(2.26b) 



to the so-called symmetric basis characterized by the 
charge vector 

^)=(l,l,... ,1)T. (2.26c) 
Example 2: The block-diagonal transformation 



V 



n—l — rn 

+1 



(2.27) 



with 1 < m < 71 < that interchanges >f„„ with >f„„, 
>c„in with — >f„„j, while it substitutes — f3„ for p„ and 
+Qni for Qn- 



III. 



EDGE THEORY 



In this Section, we study the quantum field theory for 
2N Abelian Chern-Simons fields as defined in (2.1a I or. 



equivalently, (2.2b) in a system with a boundary by fol- 



lowing a strategy pioneered in Ref s . 1491 and 1501 However, 
before relaxi ng the condition 9 $7 = 0, we decompose the 
action (2.1a) of the bulk theory into 



S :— Sp^ + Sq+ Sg, 



(3.1a) 



Sk := 1 dt 1 d^x K,^ e^^P a,,^ a,- (3.1b) 



R n 



5n := + / d< / d^x 



27r 



e'^-''a,,,9,4 , (3.1c) 



Ss:=+ j dt I d'x ^ e'^'-P a,,^ B^. (3. Id) 



Notice that we have performed a partial integration in 
Eq. (3.1c) and Eq. (3. Id) as compared to Eq. (2.1a), so 



that the gauge fields A and B enter Eq. (3.1) in an ex- 



plicitly gauge invariant form. In contrast, we are going to 
make a gauge choice for the fields with i = 1, . . . , 2N 
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to derive the gauge-invariant effective theory of the edge, 
once we have relaxed the condition 90 = 0. 

Let us choose f2 to be the upper-half plane of M-^, i.e., 



n:^{{x,y)eR^\y>0} 



(3.2) 



for notational simplicity but without loss of generality. 
Observe that under the 2N independent Abelian gauge 
transformations of the dynamical Chern-Simons fields 

a»,M ^ + d^,Xi (3.3a) 

for fjL — 0,1,2 where Xi with i ~ I,-- - , 2A^ are real- 
valued and smooth, the action S defined in Eq. (3.1) 
obeys the transformation law 



S + 5S 



with 



+ 00 



+ 00 



5S 



dt / dx (x, Ji')(t,x,0) 



(3.3b) 



(3.3c) 



boundary arising from the quasi-particle current normal 
to the boundary, i.e.. 



0= j j J^{t,X,Q). 
— oo —oo 

Observe that the stronger condition 



X,(i,x,y = 0) =0 



(3.8) 



(3.9) 



for i = 1, • • • , 2N achieves gauge invariance, i.e., SS — 0, 
without imposing condition (3.8 1. 



Now that we understand under what conditions the 



quantum field theory with the action (3.1) is gauge in 



variant with the choice (3.2) for we are ready to con 



struct the bulk-edge correspondence. To this end, we 
are going to extract from the dynamical gauge field a 
degrees of freedom that are localized on the edge dfl 
and invariant under the gauge transformations induced 
by Eqs. (3.3a I, (3.6), and (3.8) on the edge 



and 



1 



An 
e 

2^ 

s 



-K,^e"'P (9, a,-,) (t,x,y) 
Q^e^'^P {d,A^){t,x,y) (3.3d) 
, ^^S,^^"' {d,B^)it,x,y). 
The equations of motion 
K,^ e^-'P a,. , = e e'^"" d,A^ + s e^'^P (3.4) 
for the dynamical gauge field a dictate here that 



jP{t,x,y) = +^K,^eP''Pd,a^^p (3.5) 

for i = 1, • • • , 2N and /i = 0, 1, 2. Hence, the 2N compo- 
nents of the quasi-particle 3-current J7j obey the conti- 
nuity equation d^J}^ = if {d^ — d,, d^) „ — holds 



,27V and p = 0,1,2. 



"1,1] 



for any i = 1, 

We now assume that the 2iV-tuplet x is constant along 
the boundary dVl for all times, 

{d^Xi) {t, x,y^Q) = {d^Xr) {t, x,y = 0)=0 (3.6) 

for i = 1, • • • , 2N. In this case, each component Xi can 
be pulled outside the integral in Eq. (3.3c I yielding 



+ 00 +00 

6S^X^ J dt J dxj^^t,x,0). 



(3.7) 



Gauge invariance, i.e., SS = 0, is then achieved if, in 
addition to the restriction (3.6 1, we demand that there 



is no net accumulation of quasi-particle charge along the 



A. Bulk-edge correspondence 

We start by fixing the gauge of the 2N Abelian Chern- 
Simons fields through the conditions 



an = K-^Va,. 



(3.10a) 



We demand here that 1^ is a symmetric, positive definite 
2A^ X 2N matrix that satisfies 



where the 2A^ x 2N matrices 



1 



p=l,2,3 



(3.10b) 



(3.10c) 



are defined by taking the tensor product between any 
of the Pauli matrices CTj^, ctjj ^^'^ "'3 ^^^^ the unit N x 



N matrix Ijy- Condition (3.10b) guarantees that the 



gauge condition (3.10a) is consistent with reversal of time 
defined by 



aM,x,y) ^ -gP"^! a^{~t,x,y). 



(3.10d) 



Indeed, the gauge condition (3.10a) then transforms un- 
der reversal of time into 

-j:^a°i-t,x,y)=K-^VJ:,a^i~t,x,y), (3.11) 

which, upon using K^^ — —J^iK~^T,i, coincides with 
Eq. (3.10a) if and only if we impose condition (3.10b). 



Next, we use the gauge conditions (3.10a) to eliminate 
the time components Oq of the dynamica gauge fields 
from the theory. For that, observe that their equations 
of motion 



= 



San 



9x0.2 - di^i = 0, 



(3.12a) 
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which require the vanishing of their field strengths, are 
automatically satisfied if 



then follows if the 2N components $j of the vector field 
$ are smooth for i = 1, • • • , 27V. 



for 



fli = 02 = ^2$, (3.12b) 

{8^82-828^)^ = (3.12c) 



We rewrite the kinetic part (3.1b) of the action (3.1a) using the gauge conditions (3.10a| and the equations of 



motion (3.12a| and subsequently substitute the gauge fields $ defined in Eq. (3.12b) 

hoo +00 +00 



s 



K 



An 



An 



J dt J dx J dy {-alK8,a^+alV8^a,) 



— 00 


—00 





+00 


-\-oo 


+00 




1'" 


/ 


—00 


— 00 





+00 


+00 


+00 



j dt j dx j dy8^{^'^ K 8^8)^^-^'^ V 8)^8^^) 



— 00 —00 

We shall demand that <^{t^x) vanishes for \x\ oo, in which case 

+00 +00 



Sk=- 



1 

An 



^ J dt J dx ($"^ifao<9i$-$"^V^9i9i$) it,x,0) 
{8i<^f K 8o<l> - {8^^)'^ V8M {t,x,0). 



An 

—00 —00 
+00 +00 

dt J dx 

-00 —00 



(3.13) 



(3.14) 



Under the gauge transformation (|3.3a|) subject to the y = 0, i-e. 

iforms 

(3.15) 



constraints (3.6) and (3.8) the 2A^-tuplet $ transforms 
as 



<i>{t,x) $(i,cc) + x- 

The fact that x is independent of time t and space x 
implies that (a) the edge theory (3.14) is unchanged un- 
der Eq. (3.15), as anticipated, and (b) { 81^) {t, x, 0) and 
(i9o$) {t, x, 0) are unchanged under Eq. ( 3.15| ) and there- 
fore are physical degrees of freedom at the edge. Their 
dynamics are controlled by the non-universal matrix V, 
which is fixed by microscopic details of the physical sys- 
tem near the edge. 

So far, we have discussed only the kinetic part of the 
action. Let us now discuss the couplings to the external 



gauge potentials A and B given by the actions (3.1c 



and (3.1d), respectively. We assume that the external 
gauge field A is chosen so that (i) all its components are 
independent of y, i.e., 

A^{t,x,y)=A^{t,x), (3.16a) 

for /z = 0, 1, 2 and (ii) they generate the Maxwell equa- 
tions in a one-dimensional space defined by the boundary 



A2{t,x) = 



(3.16b) 



for all times t and for all positions x along the one- 
dimensional boundary y = 0. Using (i) and (ii), we can 
recast Sn as 



Sr 



dt j d^xQ^e^^Pa.^^d.A^ 



00 00 



= J dt J dx {e^''A^Q''8,^) it,x,0). 

— 00 —00 

(3.17) 

On the last line, the Levi-Civita tensor is defined for 
(1 -f 1) space and time. 
Furthermore, the very same manipulations that lead 
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to Eq. (3.17) can be carried out on Sg to deliver 

+ C30 +00 



-00 —00 



(3.18) 



Finally, the operation of time reversal stated in 
Eq. (3.10d) in the bulk reduces on the boundary to the 



transformation law 

a,it,x) = (a,$) it,x) 
T 



Eiai(-t,a;) = (a,I]i<i>) (-i,x). 



(3.19) 



The transformation law of the 2A^-tuplet $ under reversal 
of time is thus only fixed unambiguously up to an additive 
constant 2A''-tuplet. The choice 



x) ^ Si $(-i, x) + nK-^ Q, 



(3.20a) 



with 



:= o (So 



:= - (So - S3) , (3.20b) 



guarantees that at least one Kramers doublet of fermions 
exists as local fields in the edge theory, as was shown in 
Ref.Ha 



B. Fermi-Bose edge correlation functions 

Local excitations on the edge can be classified into two 
groups. There are quasiparticle excitations that carry 
rational charges and obey fractional statistics. There 
are Fermi-Bose excitations that carry integer charges and 
obey Fermi or Bose statistics. The former excitations are 
built from vertex operators of the form 



vrit,x) 



{t,x) 



(3.21a) 



that are labeled by the flavor index i = 1, - ■ ■ , 2N. The 
latter excitations are built from the vertex operators of 
the form 



(3.21b) 



that are also labeled by the flavor index i = 1, • • • , 2A''. 
Establishing the statistics under exchange obeyed by 
these vertex operators can be achieved by computing 
their correlation functions, as we now show for the Fermi- 
Bose operators. 

We shall choose for H. a disk of unit radius centered at 
the origin of the complex plane with coordinate z € C 
Thus, the boundary dfl is the unit circle centered at the 
origin of C. We are after the correlation function 



^e^ ^1^(2:1,1, %^i) X • • • X V^i^(zi % „J X • • • X V.^{z2N i, 22jv,i) ^ 

I 



(3.22) 



where the angular bracket denotes an expectation value 
using the quantum field theory with the action (3.14) 



and Q is a so-called background charge. This correlation 
function fixes the positions of particles of flavor i = 



1, • • • , 27V at the locations i, Zj . 



along the 



unit circle. We are omitting any reference to the time t 
since all Fermi-Bose vertex operators are taken at equal 
time. 

We shall use the rules that 




'w), if I = J =!,■■■ ,N, 
(^^j{z,z)^j{w,w)'^ ^ l\og(z-w), ifI = J^N + l,--- ,2N, (3.23) 

otherwise, 



where the capitalized index / = 1, • • • , 2N labels the ba- 
sis of for which the K matrix is represented by the 
diagonal matrix made of the signature of its eigenvalues 



)^ K {W~ 



(3.24) 



Observe that the linear transformation W needs neither 
be integer- valued nor have unit determinant. It is a mere 
useful device to compute the correlation function (3.22 1. 



The relationship between the co-ordinates $j and $ is 
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linear and given by 

^j = Wj^^^, /-l,---,2iV, (3.25) 
where the summation convention for the repeated small 



case indices is used. 



In order to take advantage of Eq. (3.23 1 when evaluat 



ing Eq. (3.221, we use the decomposition 



(3.26a) 



for i = 1, • • • ,27V where the matrices and T,^ were defined in Eq. (3.20b|. Under the decomposition (3.26a 



Fermi-Bose vertex operator (3.21b I occurring in the correlation function (3.22) becomes 

V^iz, -z) = e-'^.. *.(^-^") = exp (-i {Wj, eJ, $,) (z)) x exp (+i {w,, ^j) [z) 

We shall also decompose accordingly the background charge Q — + Q^. Now, 

/ . 



any 



2t ■ 2 



/ / 2W \ 
\ y 1=1 a.=l J 
I I 2N ih \ 
\ y j=l a^=l J 

11 11 \^^,<^,~\b,) [z^,a^- Z,^b,) 



X exp 



2N 



i—1 l<a. <b- <n. 



n li ( 

l<i<j<2N a. ,b,=l 



W„T.t,W, 



(-Z --Z Y" " 



(3.26b) 



(3.27) 



The role of the background charge is to guarantee "charge 
neutrality". Observe that for any pair i,j = 1, • • • , 2iV, 
the exponents 



and 



a,- 



(3.28) 
(3.29) 



can be an irrational number! Nevertheless if 2^ ^ ^ z i 
the correlation function (3.27) is single- valued, for it is 



the product of functions of the form 



/(^i,a. 1 ^7,6 . ) 



ft. 



(3.30a) 



where one verifies that 

^ij l^ij 



is integer valued. This is consistent with the fact that 

(3.30b) 
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the vertex operators (3.21b I describe either Fermions or 
Bosons. 



IV. WAVE FUNCTIONS 

The family of topological quantum field theories de- 
fined by Eqs. ( 2.2b ) and ( 2.25 ) encode the universal prop- 



Eq. (2.26l]. Second, for any representation of the univer- 



erties of a family of time-reversal symmetric fractional 
quantum liquids. 

Any connection to a microscopic realization of a time- 
reversal symmetric fractional quantum liquid whose uni- 



versal properties are captured by Eqs. (2.2b) and (2.25) 
must, however, be supplied. 

For example, the very definition of an electron op- 
erator is ambiguous for any equivalence class of topo- 
logical quantum field theories defined by Eqs. (2.2b) 



and (2.25). First, there is no unique definition of a lo- 



cal fermion in any topological quantum field theory of 
the form Eqs. (2.2b I and (2.251 that admits a hierar- 



chical representation with one type of Kramers degener- 
ate pairs of fermions, for this representation is equiva- 
lent to the symmetric representation that admits N dis- 
tinct types of Kramers degenerate pair of fermions [see 



sal data {K, Q, S) from Eq. (2.1b) that admits fermions. 



a basis set of functions must be supplied to construct a 
representation of the microscopic electron operator. This 
basis set of functions is usually provided by some refer- 
ence single-particle electron basis. 

In the context of the FQHE observed in a GaAs ac- 
cumulation layer, the basis set of functions is the single- 
particle basis of the Landau Hamiltonian describing an 
electron moving in a plane perpendicular to a uniform 
magnetic field. However, the Landau basis set and, in 
particular, the basis set for the lowest Landau level, is not 
appropriate for the recently discovered fractional quan- 
tum Hall phases in lattice models without external mag- 
netic fieldPSEH 

With this caveat in mind, we are going to con- 
struct some wave functions using the data {K, Q, S) from 
Eq. (2.1b) as the universal input and using the Landau 



wave functions spanning the lowest Landau level as the 
microscopic input. We do this out of simplicity in view 
of the elegant analytic properties of these single-particle 
functions. Hence, we choose the symmetric gauge for 
which the Slater determinant in the lowest Landau level 
is 



4- 



=1 ({^^i , Z; ) 



n 



l<A;<;<n. 

for n- electrons labeled by the flavor index i = 1, • • • ,N while it is 



fc=i 



n 



l<fc<;<n. 



k=l 



-i,fe^i,fe 
4^2 



(4.1a) 



4P 



(4.1b) 



for njY_|_| electrons labeled by the flavor index + i = 
A^-t-l,--- ,2A^. Here, 

{^i ' In. {^Ll' ' ' ' I ^i,riii ' ■ ■ ; ^i,ni } (4-2) 

denotes the complex coordinates of the particles of the 
first A^ flavors, with z denoting their complex conjugates, 
and likewise {w- , Wj }^ denotes the complex coordinates 

for the last A^ fiavors; i'^ := 0o/(27r|B|) is the square of 
the magnetic length £ in the presence of the uniform 
magnetic field of magnitude \B\, (pQ := 2 7r/e is the quan- 
tum of magnetic flux, |ri||-B| is the magnitude of the flux 
threading the disk fi, and ly- := {n-cl)^) / {\U,\\B\) repre- 
sents the filling fraction of the lowest Landau level. 

It remains to decide on the number of electron fiavors, 
a microscopic input. We shall assume that the hierarchi- 
cal (symmetric) representation corresponds to a single 
pair (A^ pairs) of microscopic fiavors of electrons form- 



ing a Kramers doublet (A^ Kramers doublets). We begin 
with the wave function for A^ = 1, in which case there is 
no distinction between the two representations. We then 
work out examples with A^ = 2 in the symmetric and the 
hierarchical representation. 



Wave function for A = 1 



We choose the universal data to be 

^+771 



K 



Q eGL(2,Z), Q=(J]e^2 



(4-3) 

for some given positive odd integer m. The spin filling 
fraction defined in Eq. (2.4) is 



1 



(4.4) 



11 



The putative ground state wave function that generahzes (see Ref. to the time-reversal symmetric case is 
the V = 1/m single-layer wave function from Laughlin 



n n 

i=l j=i+l 



X JJ^exp 

1=1 



(4.5) 



By construction, it is invariant under the operation of 
time reversal represented by 



r _ 



Wi ^ Zi, 



I,-- - ,n. 



(4.6) 



It thus realizes a time-reversal symmetric fractional in- 
compressible state. Observe that this wave function 
factorizes into an holomorphic and an antiholomorphic 
sector. Time-reversal symmetry forbids a coupling be- 
tween the holomorphic and antiholomorphic sector when 
N = 1. 



B. Wave functions in the symmetric 
representation 



We choose the universal data to be 
K 



mi n 



eGL(4,Z), Q = 
The 2x2 matrix k is given by 

e GL(2,Z) 



(4.7a) 



(4.7b) 



We impose that the integers rrii and mj are odd and posi- 
tive while the integer n is positive [n > is not restrictive 



in view of Eq. (2.27)] whereby 



m^mj - > 0, (4.7c) 
in order for k to be maximally chiral. In turn, 



A 



+d 
-d 



(4.7d) 



where the integer d > is chosen to be non-negative. 
Finally, the charge vector 



e = 



(4.7e) 



enforces the presence of 4 fermions related pairwise by 
reversal of time. We assume that the fermions with the 
charge vector Q in the topological quantum field the- 
ory represent 2N distinct flavors of electrons in a micro- 
scopic theory. For example, each flavor of electrons could 
be constrained to move with the dynamics dictated by 
the single-particle Landau Hamiltonian in its own two- 
dimensional layer in the presence of a uniform magnetic 
field pointing up for the first N layers and down for the 
next N layers. If the lowest Landau level of each layer is 
partially filled, interactions might select an incompress- 
ible ground state. The spin filling fraction defined in 
Eq. (pll) is 



2n 



d2- 



(4.8) 



The putative ground state wave function that generalizes 
the {mi,m2,n) bilayer wave function from Halperin (see 
Ref. H]) to the time-reversal symmetric case is 



({^I'^lin, ;{^2,^2}„^ ;{^«2,U)2}„J = 

*l/mi {{Zl,Zi}n, X *l/m, ({^2,22}„, \ {w2,W2}n,) 



(4.9) 



i=ii=i 

I 



W' 



2,i 



By construction, it is invariant under the operation of for i = 1,2. It thus realizes a time-reversal symmetric 
time reversal represented by fractional incompressible state. 



r _ r _ . -, 

■ ' ' W;, -J> Z;j_ , 2i = 1, • • • ,«! , 



(4.10) 
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Wave functions in the hierarchical 
representation 



We choose the universal data to be 
-A' 



K 



-A -K 



eGL(4,Z), Q 



The 2x2 matrix k is given by 

+m +1 



-1 



-P 



e GL(2,Z) 



(4.11a) 



(4.11b) 



where m is a positive odd integer and p is an even integer 
larger than zero. The 2x2 matrix A is given by 



+d 
-d 



with d any positive integer. Finally, the charge vector 



(4.11d) 



enforces the presence of 2 fermions related pairwise by 
reve rsal of time. The spin filling fraction defined in 
Eq. (pll) is 



mp - 



1 -d2- 



(4.12) 



The putative ground state wave function that generalizes 
the v = ^p_^_i single-layer wave function from Halperin 



(4.11c) (see Ref. 2) to the time-reversal symmetric case is 



n/d^./d^c, 



pn n 

i=ii=i 



(4.13) 



By construction, it is invariant under the operation of 
time reversal represented by 



r _ 



for z = 1 , • • • , p7i and 



(4.14a) 



(4.14b) 



for i — I, ■ ■ ■ , n. It thus realizes a time-reversal symmet- 
ric fractional incompressible state. 



V. SUMMARY 

In this paper, we first derived a hierarchy of FQSHEs, 
the universal properties of which are encoded by equiva- 
lence classes of BP theories of the form 



(5.1) 



The N X N invertible and integer-valued matrix >f cou- 
ples the N flavors of the dynamical gauge field 0^+^ to 



the N flavors of the dynamical gauge field a'"-*. The A^- 
tuplets a^~^^ and a^~^ also couple linearly to the external 
gauge fields A and B, respectively, through the vector 
g G Z^, where the integer g- shares the same parity as 
the integer for i = 1, • • ■ ,N. Correspondingly, there 
exists two independent conserved currents, a charge cur- 
rent associated to the gauge field a^"''-' and a spin current 
associated to the gauge field a^~\ 

Time-reversal symmetry implies the vanishing of the 
charge Hall conductivity 



Ztt 



The non- vanishing spin filling fraction 



0^ X ^ Q 



can be interpreted as the spin Hall conductance 



ZTT 



(5.2) 



(5.3a) 



(5.3b) 



if the U(l) conservation law associated to the cur- 
rent of a^~^ arises microscopically from a residual spin- 
1/2 U(l) (easy plane XY) symmetry. The topological 
ground state degeneracy, if two-dimensional space fl has 
a toroidal geometry. 



(det k)"^ 



(5.4) 
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is always the square of an integer as a consequence 
of time-reversal symmetry. Equivalent pairs (>£■, g) and 
g'), as defined by Eq. ( 2.25[ ), share the same spin 
Hall conductivity and topological degeneracy. 

The theory (5.1) is topological when two-dimensional 
space r2 has no boundary, i.e., the Hamiltonian den- 



sity associated to the Lagrangian density (5.1) vanishes. 
This is not true anymore if the boundary is a one- 
dimensional manifold. We have shown that imposing 
gauge invariance delivers a gapless theory with all ex- 
citations propagating along the boundary d^l. These ex- 
citations can all be constructed out of N pairs of counter- 
propagating chiral bosons whose non-universal velocities 
along the boundary dO, derive from a gauge-fixing con- 
dition. The stability of this edge theory to the (time- 
reversal symmetric) breaking of translation invariance 
along the boundary (including the breaking of the spin 
conservation law associated to the spin vector S) was 



studied in Ref. The correlation functions for the 

Fermi-Bose excitations along the edge were computed 
and shown to be a product over the functions (|3.30a|). 



Finally, we have proposed a time-reversal symmetric 
counterpart to the hierarchy of wave functions that have 
been proposed in the context of the FQHE by way of 
few examples, the — 1/m, = pj (mp -I- 1 — (P), and 



(mi 



2n) / {mim2 — n + d ) sequences. 



Acknowledgments 



We gratefully acknowledge useful discussions with 
Maurizio Storni. This work was supported in part by 
DOE Grant DEFG02-06ER46316 and by the Swiss Na- 
tional Science Foundation. 



^ R. B. Laughlin, Phys. Rev. Lett. 50, 1395 (1983). 

^ B. I. Halperin, Helv. Phys. Acta 56, 75 (1983). 

^ F. D. M. Haldane, Phys. Rev. Lett. 51, 605 (1983). 

* B. Halperin, Phys. Rev. Lett. 52, 1583 (1984). 

^ R. B. Laughlin, Surf. Sci. 142, 163 (1984); ibid. 142, A239 
(1984). 

" S. M. Girvin, Phys. Rev. B 29, 6012 (1984). 

A. H. MacDonald and D. B. Murray, Phys. Rev. B 32, 
2707 (1985); A. H. MacDonald, G. C. Aers, and M. W. C. 
Dharma-wardana, Phys. Rev. B 31, 5529 (1985). 

* J. K. Jain, Phys. Rev. Lett. 63, 199 (1989); Phys. Rev. B 
40, 8079 (1989); ihid. 41, 7653 (1990). 

^ N. Read, Phys. Rev. Lett. 65, 1502 (1990). 
^° X. G. Wen and A. Zee, Nucl. Phys. B 15, 135 (1990). 
" B. Blok and X. G. Wen, Phys. Rev. B 42, 8133 (1990); 

ibid. Phys. Rev. B 43, 8337 (1991). 
^2 J. Frohlich and A. Zee, 1991, Nucl. Phys. B 364, 517 

(1991) . 

X. G. Wen and A. Zee, Phys. Rev. B 46, 2290 (1992) 
" Xiao-Gang Wen and A. Zee, Phys. Rev. Lett. 69, 1811 

(1992) . 

J. FrohUch and E. Thiran, J. Stat. Phys. 76, 209 (1994). 
R. B. Laughlin Phys. Rev. B 23, 5632 (1981) 

B. I. Halperin, Phys. Rev. B 25, 2185 (1982). 
X.-G. Wen, Phys. Rev. Lett. 64, 2206 (1990). 
X.-G. Wen, Phys. Rev. B 43, 11025 (1991). 

2° X.-G. Wen, Phys. Rev. B 44, 5708 (1991). 

21 J. Frohlich and T. Kerler, Nucl. Phys. B 354, 369 (1991). 

22 X. G. Wen, Phys. Rev. B 40, 7387 (1989). 

23 X.-G. Wen, Int. J. Mod. Phys. B4, 239 (1990). 

2'* M.C. Diamaritini, P. Sodano, C.A. Trugenberger, Eur. 

Phys. J. B53, 19 (2006). 
2^^ N. Read and S. Sachdev, Phys. Rev. Lett. 66, 1773 (1991). 
2^^ X. G. Wen, Phys. Rev. B 44, 2664 (1991). 
2^ Christopher Mudry and Eduardo Fradkin, Phys. Rev. B 

50, 11409 (1994). 
2* T. Senthil and Matthew P. Fisher Phys. Rev. B 62, 7850 

(2000). 

2^ R. Moessner, S. L. Sondhi, and Eduardo Fradkin, Phys. 
Rev. B 65, 024504 (2001). 



''^ M. Freedman, C. Nayak, K. Shtengel, K. Walker, and Z. 
Wang, Ann. of Phys. (N.Y.) 310, 428 (2004). 

31 T. H. Hansson, V. Oganesyan, and S. L. Sondhi, Ann. 
Phys. 313, 497 (2004). 

32 S.-P. Kou, X.-L. Qi, and Z.-Y. Weng, Phys. Rev. B 71, 
235102 (2005). 

33 S.-P. Kou, M. Levin, and X.-G. Wen, Phys. Rev. B 78, 
155134 (2008). 

3'' C. Xu and S. Sachdev, Phys. Rev. B 79, 064405 (2009). 
35 C. L. Kane and E. J. Mele, Phys. Rev. Lett. 95, 226801 
(2005). 

3** C. L. Kane and E. J. Mele, Phys. Rev. Lett. 95, 146802 
(2005). 

3'^ B. A. Bernevig and S.-C. Zhang, Phys. Rev. Lett. 96, 
106802 (2006). 

3** B. A. Bernevig, T. L. Hughes, and S.-C. Zhang, Science 

314, 1757 (2006). 
3^ M. Konig, S. Wiedmann, C. Briine, A. Roth, H. Buhmann, 

L. W. Molenkamp, X.-L. Qi, and S.-C. Zhang, Science 318, 

766 (2007). 

M. Levin and A. Stern, Phys. Rev. Lett. 103, 196803 
(2009). 

G. Y. Cho and J. E. Moore, Ann. Phys. 326, 1515 (2011). 
■12 G. Y. Cho and J. E. Moore, Phys. Rev. B 84, 165101 
(2011). 

'13 X.-L. Qi, Phys. Rev. Lett. 107, 126803 (2011). 

T. Neupert, L. Santos, S. Ryu, C. Chamon, and C. Mudry, 
Phys. Rev. B 84, 165107 (2011). 

*5 M. O. Goerbig, arXiv:1107.1986 (unpublished). 

*® This terminology has no deep meaning. It arose historically 
in the context of topological field theories and refers to 
the symbols h and / used to represent a p-form 6 and the 
curvature / of a (d — p)-form in [d + 1) -dimensional space 
and time.'^ 

^'^ M. Blau and G. Thompson, Ann. Phys. 205, 130 (1991); D. 
Birmingham, M. Blau, M. Rakowski, and G. Thompson, 
Phys. Rep. 209, 129 (1991). 
A. Blasi et at., e-print arXiv: 1106.4641. 
S. Elitzur, G. Moore, A. Schwimmer, and N. Seiberg, Nucl. 
Phys. B 326, 108 (1989). 



14 



®° Xiao-Gang Wen, Advances In Physics 44, 405 (1995). 
®^ T. Neupert, L. Santos, C. Chamon, and C. Mudry, Phys. 
Rev. Lett. 106, 236804 (2011). 

D. N. Sheng, Z. Gu, K. Sun, and L. Sheng, Nature Com- 
munications 2, 389 (2011). 



Y.-F. Wang, Z.-C. Gu, C.-D. Gong, D. N. Sheng, Phys. 
Rev. Lett. 107, 146803 (2011). 

N. Regnault and B. A. Bernevig, e-print arXiv:1105.4867. 



